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ABSTRACT 

o : 

The water maser lines in the Seyfert nucleus of the Circinus galaxy vary on time 
scales as short as a few minutes. The amplitude of one line more than doubled in ~ 10 
minutes, reaching ~ 37 Jy, which corresponded to an increase of ~ 6 L , assuming 
isotropic emission, in a single maser feature on a size scale of about 1 AU, based on 
f^, ' light-travel time. Other lines vary by up to about 30% on similar time scales. The 

variability is at least two orders of magnitude more rapid than any observed for other 
— ■ Galactic or extragalactic water masers. The intensity changes cannot be attributed 

easily to a mechanism of intrinsic fluctuations. The variability may be the result of 
, strong interstellar diffractive scintillation along the line of sight within our Galaxy. 

(— | ! This would be the first example of diffractive scintillation for any source at 22 GHz 

Q-f and for any source other than a pulsar. However, only the very shortest timescales 

for interstellar scintillation, obtained from pulsar observations and scaled to 22 GHz, 

■ correspond to the observed maser variability. Alternatively, the intensity changes may 
be a reaction to fluctuations in compact background or radiative pump sources and 

■ thereby may be related to variability of the central engine. 

The maser spectral features symmetrically bracket the systemic velocity of the 
galaxy, with components red- and blue-shifted by about ±(100 — 200) km s . The 
spectrum of the Circinus maser is similar in some respects to that of the maser in 
NGC 4258, which probably traces a molecular disk rotating around a supermassive 
object. VLBI observations could reveal whether the maser source in the heart of the 
Circinus galaxy is part of a similar dynamical system. 
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1. Introduction 



Renewed interest in extragalactic water masers has been stimulated by observations of the 
maser m the nucleus of NGC 4258 QMiyoshi et al. 1995Q . Maser emission is present at the systemic 
velocity of the galaxy and in satellite lines offset by about ±900 km s _1 . The emission arises from 
a thin, nearly edge-on molecular disk that is in Keplerian rotation around a massive central object 
(see also jMoran et al. 1995| ). These studies have established, to high accuracy, the binding mass 
and the geometry of this system on sub-parsec scales. 

The Circinus galaxy (CG) is a Seyfert 2 galaxy at a distance of about 4 Mpc that has long 
been known to harbor a water maser in its nucleus ( |Gardner fc Whiteoak 1982 ). The nucleus 



shows a one-sided ionization cone ( Marconi et al. 1994 ) and bipolar radio lobes that extend about 



200 pc along the minor axis of the galaxy ( Elmouttie et al. 1995 ). Excitation of coronal lines 
within the inner 10 pc is consistent with a photoionizing central continuum source, rather than 
with a starburst, though there are signs of an old starburst within the innermost 4 pc ( pliva et 
al. 1994). A heavily obscured hard (2-10 keV) X-ray source ( Matt et al. 1996| ) and a compact 
radio synchrotron source also lie in the nucleus ( pNforris fc Forbes 1996 ). 



The water maser spectrum of the CG is similar to that of NGC 4258 in two ways. First, 
there are two complexes of narrow lines that bracket the systemic velocity of the galaxy. 
Second, there is a distinct complex of emission lines close to the systemic heliocentric velocity 
of 416 km s _1 ( Lehnert fc Heckman 1995 ). (We use the radio definition of radial velocity 
throughout.) The red-shifted lines in the CG were discovered by Gardner & Whiteoak (1982), 
but the blue-shifted lines have only recently been detected by us (reported in this paper), and by 
Nakai et al. (1995) and Braatz, Wilson, & Henkel (1996), independently. 

Whiteoak & Gardner (1986) observed that the shape of the CG maser spectrum varies 
substantially in about one month. The relative intensities of spectral components in the NGC 4258 
maser can also vary by tens of percent in a few weeks ( prccnhill et al. 1995 ). The maser in 
IC 10 can vary by 60% in less than one day (Argon et al. 1994), which may be a consequence 
of interstellar scattering. We note that the most rapid known flares in Galactic interstellar H2O 
masers have time scales of one to ten days (e.g., Liljestrom 1993; Pacheco fc Scalise 1993| ; paschick^ 
Burke, & Spencer 1977). In this Letter, we report the measurement of a finely sampled time-series 
of spectra, which has revealed variability on time scales at least two orders of magnitude shorter 
than those known for any other water maser. We also analyze the morphology of these spectra in 
the context of an edge-on disk geometry. 
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2. Observations 

The observations were carried out at the Australia Telescope National Facility's Parkes radio 
telescope in 1995 October, as part of a search for new water maser sources in 28 galaxies (z < 0.1) 
with high X-ray absorption columns or excess infrared emission accompanied by compact radio 
continuum emission. We used a newly commissioned 22 GHz receiver, with two orthogonal linear 
polarization channels. From a tipping scan, we estimate that the receiver temperature was 95 K, 
and that the atmospheric opacity at zenith was 0.06. The telescope sensitivity was about 5.7 Jy 
K _1 , based on a flux density of 21.2 Jy for Virgo A. Uncertainties in pointing are responsible for 
up to 10% loss in gain. We adopted the gain curve of Bourke (1994), increased by about 10% at 
low elevation angles as indicated by our observations of Virgo A. We observed in position-switching 
mode, with four 16 MHz bandpasses tuned to cover the range of emission from 195-735 km s _1 . 
The 4096-channel digital autocorrelator provided a channel spacing of 15.625 kHz, or 0.21 km s _1 . 



3. The observed spectrum and variability 

The spectrum of the CG water maser shows two complexes of strong lines between 510 and 
610 km s _1 , and between 240 and 310 km s _1 (see Fig. 1). We refer to these as the red and the 
blue satellite lines. Spectra taken in previous years show similar velocity ranges for the complexes 
of satellite lines ( [Whiteoak fc Gardner 1986| ; |Nakai et al. 1995| ). Fi gure 2 shows the relative flux 



density variations for the three strongest maser lines on time scales as short as a few minutes. The 
567.25 km s _1 feature more than doubled in strength to about 37 Jy in only 10 minutes between 
15 and 16 hours LMST on day 289. On the assumption of isotropic emission this represents an 
increase from 4.8 to 11 L Q . At other times the strength of this feature varies by about 50%, while 
the other two features vary upto about 30%, on similar time scales. We conservatively adopt a 
characteristic time scale for variability, t (, s , of ~ 300 s. 

We feel certain that instrumental effects are not responsible for the apparent maser 
fluctuations. The variations of different spectral features are not correlated. Moreover, 
observations of the CG and NGC 4945 masers were made within hours of each other, and no 
variability was detected in NGC 4945. Apparent variability might be a consequence of telescope 
pointing errors if the maser features were distributed over an area comparable to the 1(3 primary 
beam. However, spectra taken at offset positions show that the region of maser emission is 
relatively compact, < 0'3, as expected for an extragalactic nuclear maser. False variability might 
also arise from under-sampling of the spectral lines but in our case the half-power widths (FWHM) 
of the maser features are at least four spectral channels. 

The FWHM of the narrowest satellite spectral line is about 0.9 km s . The satellite lines 



seen in NGC 4258 are about 1.0 km s 1 wide (Miyoshi et al. 1995), and the narrowest water maser 



lines ever observed are about 0.3 km s 1 in IC 10 (Argon et al. 1994). Emission from masers in 



the presence of a magnetic field may be linearly polarized (Watson 1994a). We used our dual 
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Fig. 1. — Time-averaged spectrum of the water maser emission from the nucleus of the Circinus 
galaxy. The heliocentric systemic velocity of the galaxy, 416 km s _1 , is indicated by the arrow. The 
dashed lines indicate the divisions between the four receiver bandpasses. The on-source integration 
time for the two bandpasses that contain the two strong complexes of lines was 350 minutes (la 
noise = 0.05 Jy). The integration time for the other bands was 100 minutes (la noise = 0.1 
Jy). The spectrum was Hanning smoothed and has a channel spacing of 0.21 km s _1 . The three 
spectral features examined in Figure 2 are indicated by the bars along the top axis, (insert) 
Spectrum obtained separately that shows more clearly the emission near the systemic velocity (as 
a consequence of better baseline removal). The on-source integration time was 60 minutes and the 
channel spacing was 0.8 kms- 1 . The gain calibration is accurate to ~ 10 % and includes corrections 
for atmospheric attenuation, antenna gain variations, and average pointing error. 



polarization observations of the satellite lines to construct Stokes Q spectra of the source. There 
is no statistically significant linear polarization at a level of 0.5 Jy (4% for a 12 Jy line). 



4. Origins of Variability 

The maser variability is extraordinary because a light travel time of a few minutes implies a 
characteristic scale of about 1 AU. Masers may be visualized typically as filamentary structures, 
viewed end-on, which mark local enhancements in density, velocity coherence, or excitation 
conditions flElitzur 1992 ). Large fluctuations in emission rate are caused by changes in the maser 



medium over lengths greater than the maser gain length, £, which is the characteristic length 
over which the line-center opacity changes by ~ 1. For typical conditions, I S> 1 AU (e.g., a 
maser pump efficiency of 1%, a fractional water abundance of 10 -4 , and a temperature of 400 K 



as in Reid fc Moran 1988|). To be responsible for the observed maser variability, changes in the 



maser medium must propagate at about the speed of light, and we know of no suitable driving 
mechanism. 
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Fig. 2. — Radio light curves for three of the strongest spectral components, (top) 567.25, (middle) 
562.00, and (bottom) 299.91 km s~ x . Instrumental effects cannot be responsible for the time 
variability. Before 16 LMST data were taken in 5 minute integrations, and after 16 LMST in 
1 minute integrations. Data gaps occurred when the telescope was off source. The error bars reflect 
the formal uncertainties from fitting Gaussian line profiles to the spectra. It is possible that an 
error in the calibration at the lowest elevations causes a slow upward trend in flux densities after 
19 LMST, on the order of 10%. 



There are three mechanisms by which the maser may vary in response to processes external 
to the maser medium: interstellar scintillation, amplification of a time-variable background source, 
and excitation by a time-variable radiative pump source. The last two mechanisms assume the 
presence of a compact (1 AU) source in the vicinity of the maser. Large amplitude fluctuations on 



time scales close to 10 s have been observed towards nuclei at X-ray energies ( Mushotzky, Done 



fc Pounds 19~93[ ), as in NGC 4051 flGuainazzi et al. 1996|) and Mrk 766 (|Leighly et al. 1996; ). 

Interstellar scintillation (ISS) is caused by scattering in an ionized component of the Galactic 
interstellar medium, which may be important because the CG lies at a Galactic latitude of only 
—3.8° (longitude 311°). We consider scintillation caused by a thin, turbulent scattering screen at 
distance D that is moving with respect to the observer at velocity v. The turbulence is usually 
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described by a Kolmogorov power law ( Armstrong, Cordes, fc Rickett 1981| ) and the scintillation 



is divided into two regimes, weak and strong. For a point source, the time scale of weak ISS 
is approximately the crossing time of the Fresnel scale, £p ~ \JD\/2it, where A is wavelength 
([Narayan 19931 ). For D = 1 kpc, A = 1.35 cm, and v = 100 km s -1 , £ F ~ 2.5 x 10 10 cm and 



tp = £ F /v rj 4000 s ( |Rickett 1990| ), about an order of magnitude greater than t Q b s . For strong ISS 



the contributing turbulent length scales obey l sc in < £f- (Strong scintillation may be diffractive or 
refractive. We do not consider the second because the relevant length scales are longer.) Diffractive 
scintillation at 1 GHz with time scales of t SC i n > 20 s is observed toward pulsars for which D > 1 
kpc and 2° < b < 5° flCordes 1986|) . For Kolmo gorov turbulence, t sc i n oc v , where v is frequency. 



Goodman & Narayan (1985) show that for non-Kolmogorov scattering screens t sc i n oc v ■ . In 
either case, only the fastest observed scintillation at 1 GHz, t sc i n ~ 10 s, corresponds to sufficiently 
rapid variations at 22 GHz. However, the true angular size of a scintillating maser feature, 9 m , 
must also be less than 9 sc i n = £ sc i n /D because otherwise the scintillation is slowed by a factor of 
(Gm/dsdn) and the amplitude is similarly decreased ( |Narayan 1993| ). Hence, for strong scintillation, 
8m < Gscin < Of, which corresponds to < 7 AU at a distance of 4 Mpc, where 6f = If/D. The 
scintillation time scale for a scattering screen in the AGN, within < 1 kpc of the maser, may be 
< tobs , but the amplitude of the variations will be substantially reduced because the angular size 
of the maser source, viewed from the screen, will be large. 

The time scale of variability for a maser that amplifies a background source is not limited 
by the extent of the maser. The maser medium forms an image of the background source, and 
the angular sizes of the two are the same (see Haschick et al. 1990| ). For the CG, the background 



source must be less than a few light-minutes in size, and the apparent maser cross-section would 
be about 0.25 /ias or 1 AU. We note that substantial variability may be possible only if a fraction 
of the maser is unsaturated since the emission rate of a completely saturated maser is fixed by the 
pump rate and beam angle. In this case, the beam angle is set by the solid angle of the background 
source as seen from the maser (see |Haschick et al. 1990 ). 



The variability time scale of a radiatively pumped maser is also independent of its physical 
extent, apart from optical-depth effects. The pump luminosity, which may be supplied directly 
by the central engine, must be large enough that the pump rate exceeds the maser emission rate. 
For Galactic interstellar water masers, radiative pumping by stellar radiation generally cannot 
support observed emission rates however, substantially higher luminosities are present in AGN. 
For a single pump transition (at wavelength X p ), F up > F urn (Qf,/Q m )(x m /x p ), where F up and F l 
are the pump and maser flux densities, respectively, Qf, is the maser beam solid angle, Q m is the 
solid angle of the maser viewed from the pump, and x p and x m are the fractional bandwidths 



//of 



of the pump and maser, which are similar ( Reid fc Moran 1988] ). For N simultaneous pump 



transitions the requisite F vp is reduced proportionately. We assume the central engine has a flat 
spectral energy distribution {vFy) between ultraviolet and X-ray energies and estimate F vp at \ p 
from the inferred intrinsic X-ray luminosity of the central engine. Matt et al. (1996) obtain a 
luminosity of 4 x 10 41 / -1 erg s _1 from the observed X-ray spectrum, assuming that a fraction, 
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/ = (Oij/27r), is reflected by a portion of the obscuring circumnuclear torus into the line of 
sight, where £1r is the solid angle of the reflecting material seen from the central engine. Since 
the infrared luminosity of the nucleus is ~ 5 x 10 43 erg s _1 ( [Moorwood &: Glass 1984 ), / is on 



the order of a few percent. Based on order of magnitude estimates for the parameters involved, 
radiative pumping requires that X p ~ 0.6A r ~ 1 (^-)(^j^)(o^-)( 4 ^ pc ) 2 , where X p is expressed in 
units of microns, and ~ 1 ( [Reid fc Moran 1988| ). Hence, a radiative maser pump may be 

feasible if it consists of on the order of 10 transitions at ultraviolet wavelengths. 

Several other specific mechanisms for maser variability have been proposed in the literature, 
though none of them explains the observed time scale. First, in partially saturated masers, a 
radiative instability can cause variations with a period of L/c, where L is the maser filament 
length ( |5cappaticci &: Watson 1992 ). However, strong maser emission requires many gain lengths 



or L/c ^> t bs- Second, masers may amplify each other, and an outburst can occur when one maser 



drifts in front of another with the same line-of-sight velocity ( Deguchi fc Watson 1989 ). Third, 



rotation of a maser beam across the line of sight can mimic rapid intensity variations ( Norris & 



Booth 1981). However, for the latter two mechanisms, maser lengths of 3> 1 AU and transverse 



velocities of < 100 km s 1 require maser cross-sections of < 10 10 cm, which are improbably small. 



5. Structure of the maser source 

Based on the morphology of the H2O maser spectrum, we propose that the maser emission 
arises from a nearly edge-on Keplerian disk, as in NGC 4258. The satellite lines arise on a disk 
diameter perpendicular to the line of sight, and the systemic emission arises on the near side of 
the disk, close to the line of sight to the dynamical center. These locations naturally possess the 
greatest line-of-sight velocity coherence. For Keplerian rotation, the outer radius of the disk is 
about four times the inner radius because the velocities of the satellite lines at these radii, relative 
to the systemic velocity, differ by a factor of two (i.e., 100 and 200 km s _1 , respectively). The 
binding mass for the disk in solar units is 10 6 ro.i sin 2 i, where ro.i is the disk inner radius in 
units of 0.1 pc, and i is the inclination of the disk rotation axis. A central engine of < 10 7 Mq, 
with a gravitational radius of < 3 x 10 12 cm, can vary on time scales of 100 s (cf., Mushotzky et 
al. (1993)), as may be required by the observed maser variability. We note that for rotation speeds 
of 200 km s _1 and ro.i > 0.04, the line-of-sight centripetal acceleration of systemic maser features 
in the spectrum will be < 1 km s _1 yr _1 (cf., |Greenhill et al. 1995| ). 



6. Conclusions 

We have observed extremely rapid variability in the spectrum of the nuclear water maser in 
the Circinus galaxy. In principle, fluctuations on time scales of a few minutes may be explained 
by ISS in our Galaxy, variability of background emission that is being amplified, or variability 
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of a radiative pump. While only the shortest measured scintillation time scales for pulsars, 
extrapolated to 22 GHz, are acceptable, the sensitivity limits of the pulsar observations may hide 
still more rapid diffractive ISS. In addition, the scattering medium at low latitude is clumpy, 
which may also allow more rapid scintillation along the particular line of sight to the CG. If the 
maser fluctuations are caused by diffractive ISS, then the maser is the first example for any source 
at 22 GHz and for any source other than a pulsar. We note that the only other extragalactic 
water maser source reported to display significant variations over < 1 day (IC 10) also lies at low 
latitude, 3° flArgon et al. 1994| ). 

Since the maser spectrum consists of independently varying components, the alternative 
background amplification and radiative pump schemes require either several compact fluctuating 
sources or variable pump illumination that depends on location or Doppler velocity. Nonetheless, 
if the CG maser directly reflects the intrinsic variability of the central engine on time scales of 
a few minutes, then the maser will be a powerful probe of this heavily obscured AGN. Rapid 
variability has been observed exclusively in X-ray wavebands for other AGN but this does not 
necessarily preclude the possibility that emission from central engines also varies about as rapidly 
at other wavelengths. 

Scintillation and background amplification models place stringent upper limits of 1 - 10 AU 
on the apparent sizes of the maser components, which implies very large brightness temperatures, 
> 10 16 K. The general symmetry of the maser spectrum with respect to the systemic velocity is 
similar to that of NGC 4258, in which the masers lie in a Keplerian molecular disk bound by a 
massive central engine. In the absence of VLBI data, we propose that the same model may be 
applicable to the CG maser. The highly ordered velocity field of a thin Keplerian disk naturally 
supports highly elongated masers, which Watson (1994b) argues is a necessary condition to permit 
the inferred high brightness temperatures 

We thank D. Backer, C. Gwinn, R. Manchester, R. Narayan, and W. Watson, for helpful 
discussions. We are also grateful to Euan Troup, Harry Fagg, and Warwick Wilson for their special 
assistance at Parkes. 
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